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Research performed during the past decade revealed an important role of symmetry energy in the equation 
of state (EOS) of strange quark matter (SQM). By introducing an isospin-dependent term into the quark mass 
scaling, the SQM stability window in the equivparticle model was studied. The results show that a sufficiently 
strong isospin dependence Cj can significantly widen the SQM region of absolute stability, yielding results that 
simultaneously satisfy the constraints of the astrophysical observations of PSR J1614-2230 with 1.928 + 0.017 
Mo and tidal deformability 70 < A,.4 < 580 measured in the event GW170817. With increasing Cy, the 
difference between the u, d, and s quark fractions for the SQM in 6—equilibrium becomes inconspicuous for 
C > 0, leading to small isospin asymmetry 6, and further resulting in similar EOS and structures of strange 
quark stars (SQSs). Moreover, unlike the behavior of the maximum mass of u-d QSs, which varies with C'i 
depending on the sign of the parameter C, the maximum mass of the SQSs decreases monotonously with 


increasing C1. 
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I. INTRODUCTION 


At sufficiently high baryon number densities, strange quark 
matter (SQM) [1-4] with deconfined u, d, and s quarks is ex- 
pected to be the true ground state of quantum chromodynam- 
ics (QCD). Owing to its possible implications in many fields, 
e.g., nuclear physics, astrophysics, and cosmology, it has gar- 
nered significant attention in the last few decades. Theoret- 
ically, SQM is a system dominated mainly by strong inter- 
actions, which can be studied in detail using QCD. Unfor- 
tunately, owing to the well-known quark confinement phe- 
nomenon at relatively low Densities, which could be relevant 
to neutron-star interiors and quark stars, perturbative QCD 
calculations of SQM are not reliable [5]. In addition, al- 
though the ab initio lattice QCD calculations were very suc- 
cessful in the finite-temperature regime with a zero chemical 
potential (baryon number density), they still suffer from the 
famous intractable Fermion sign problem [6]. Therefore, re- 
searchers typically resort to phenomenological models, such 
as the density-dependent model [7—11], quark cluster model 
[12, 13], quasi-particle model [14-18], Nambu-Jona-Lasinio 
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model [19], and so on. In Refs. [20, 21], strong magnetic 
field effect on the properties of SQM had been investigated, 
and numerous interesting results have been obtained in other 
studies [22-25]. 

To study the SQM properties using phenomenological 
models, one of the most important issues is determining rea- 
sonable values of the model parameters. Given the conjec- 
ture that SQM is absolutely stable, the most commonly used 
constraint is that, for SQM, the energy per baryon should be 
smaller than 930 MeV, while for the bulk u-d quark matter 
(u-d QM) it should be larger than 930 MeV. Moreover, sig- 
nificant progress has been made in the observations of global 
properties of compact objects, such as mass, radius, and tidal 
deformability. For example, the measurements of PSR J1614- 
2230 and PSR J0348+0432 several years ago precisely deter- 
mined masses of 1.97 + 0.04 Mo [26] (a more precise mass 
measurement yielded 1.928 + 0.017 Mo [27] with a 68% 
credible level) and 2.01 + 0.04 Mọ [28], respectively. In 
addition, the more recently measured gravitational mass of 
PSR JO740 + 6620 reached 2.14+0:79 Mo with 95.4% cred- 
ibility, likely making it the the largest compact star observed 
up to date [29] (the updated result for PSR JO740 + 6620 is 
2.081032 Mo [30] with 68.3% credibility). Therefore, these 
observations of pulsar masses have been widely used to con- 
strain the EOS of compact stars and the onset of the phase 
transition from hadronic matter to quark matter in hybrid stars 
[31, 32]. 

In addition, with the coming of multi-messenger observa- 
tions of neutron stars, owing to the detection of gravitational 
waves for events GW170817 [33] and GW190814 [34], addi- 
tional stringent constraints on the EOS of compact star mat- 
ter have been appeared [35-37]. Specifically, it was deter- 
mined that for GW190814, the newly discovered compact 


binary merger possesses a secondary component with mass 
as large as 2.50—2.67 Mo at the 90% credibility level. If it 
is confirmed to be a compact star, rather than a light black 
hole, it will be possible to rule out many EOSs. Moreover, 
this detection of gravitational waves for GW170817 offers 
a novel opportunity to directly probe the properties of mat- 
ter under extreme conditions found in the interior of com- 
pact stars, yielding improved estimates of tidal deformability 
Aya = 1904430 [38]. This measurement of tidal deforma- 
bility has been widely used for constraining the properties of 
compact star matter [39—46]. 


Recently, an interesting proposition was made, suggest- 
ing that, instead of the SQM, the u-d QM may be the true 
ground state of QCD [47]. According to this proposition, it 
has been found that the properties of the u-d QS are in ac- 
cordance with various astrophysical observations [48], where 
quark symmetry energy effects are included in the EOS of 
the u-d QM. In practice, it has been found that symmetry en- 
ergy plays an important role in determining the properties of 
isospin asymmetry matter [49-53]. For nuclear matter, al- 
though the symmetry energy behaviors have been well deter- 
mined below and around the saturation density, it remains a 
big challenge to determine the density behaviors of the sym- 
metry energy at suprasaturation densities, where phase tran- 
sitions of strong interaction matter may occur, and various 
new forms of matter may appear. In fact, due to the different 
symmetry terms in different phases, isospin effects are rather 
significant [54]. In this study, considering the isospin effects 
of the SQM, we investigated the stability window and sym- 
metry energy of SQM, and employed some of the constraints 
mentioned above in the equivparticle model [23] to examine 
whether the SQS properties are consistent with the reported 
astrophysical observations. Itis found that inter-quark isospin 
effects are necessary for reconciling the SQSs obtained in the 
equivparticle model with the aforementioned constraints. 

The remainder of this paper is organized as follows. Sec- 
tion ILA briefly describes the equivparticle model, where 
quark mass scaling with isospin effects is introduced. Then in 
Sect. II B, the equations for calculating the tidal deformabil- 
ity of compact stars are presented. In Sect. III, the numerical 
results of the properties of SQM and the structures of SQSs 
are given and discussed. Finally, we conclude this paper with 
a brief summary in Sect. IV. 


Il. EQUATIONS FOR SYMMETRY ENERGY AND TIDAL 
DEFORMABILITY 
A. Equivparticle model with isospin asymmetry 


In the equivparticle model, by introducing an effective 
quark chemical potential už, the thermodynamic potential 
density of the quark matter system can be written as [23] 
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where f(x) = xx? —1(2? — 3) + žarch(zx), gi(=6 for 
quarks and 2 for electrons) is the particle degeneracy, m; 
is the isospin density-dependent equivalent quark mass, and 
vi = \/ už? — m? is the Fermi momentum of particle species 
i. Using the basic thermodynamic relation p; = —dQo/duž, 
the particle number density p; is readily obtained 
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Regarding the equivalent quark mass m;, different types 
of quark mass scaling were constructed to account for strong 
inter-quark interactions [7, 10, 23, 55]. Conventionally, in 
the equivparticle model, m; is a function of the baryon num- 
ber density æ. To include the isospin effect on the EOS 
of quark matter, the isospin asymmetry term ô = 3(pa — 
Pu)/(Pa + Pu) is introduced into the mass scaling frame- 
work. In Ref. [56], m; is parameterized as mj; (pp,6) = 
mio + Dp” — T,0Dipte with 7; being the third com- 
ponent of the isospin for quark flavor 7. In this study, we used 
the parameterization that was recently proposed in Ref. [48], 
i.e., 


Mil po, Ò) = Mio + Mii ( P, 6) 
= Mio + Dg + Ca + Cyd" pp, (2) 


where m;o is the quark current mass, D is the confinement 
parameter indicating the linear confinement effect, C repre- 
sents the one-gluon exchange interaction for C < 0 and the 
leading-order perturbative interaction for C > 0, and parame- 
ter Ch; (Ciu = Cla, Cis = 0) is used for adjusting the strength 
of the isospin-dependent equivalent quark mass and thus the 
magnitude of the quark-matter symmetry energy. In addition, 
this equivalent quark mass satisfies the exchange symmetry 
between u and d quarks in quark matter in the absence of 
inter-quark electromagnetic interactions. 

The energy density and pressure of the SQM are, respec- 
tively, given as 
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where g(x) = xyz? — 1(2x? — 1) — arch(x) and 
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Based on the definition of the isospin asymmetry ô, the 
derivatives of ô with respect to pu and pg are, respectively 
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Accordingly, the summation of the index 2 in the third term 
on the right-hand side of Eq. (4) is zero, which leads to the 
following simple pressure form: 
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To obtain the EOS of the SQM, one has to solve the fol- 
lowing equations, which are respectively the electric-charge 
neutrality conditions 


2Pu — Pd — Ps — 3Pe = 0, (9) 


the conservation of the baryon number density 


1 
Po = 3(Pu + pa + ps), (10) 
and the 6-equilibrium condition 


Hu + He = Md = Hs- (11) 


It should be noted that the chemical potentials in Eq. (11) 
are real ones, instead of the effective ones. The relationship 
between the real chemical potential 4; and effective chemical 
potential už in the equivparticle model is 
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Using Eq. (12), Eq. (11) can be rewritten as 
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Using Eqs. (9) and (10), Eqns. (14) and (15) can be con- 
sidered as a set of equations for p,, and pg. By solving these 
equations for a given pp, we obtain pu and pa, and then, ac- 
cordingly, the EOS of the SQM. 

There are two methods for obtaining the symmetry en- 
ergy. The first approach uses the parabolic law, approximately 
given by [56, 57] 
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Esym(Pv, Ps) > g E (Po ô= 3, ps) = E(P, é= 0, ps)]; 
(16) 


which implies that the symmetry energy is an estimate of the 
energy cost required for converting all u quarks in the sym- 
metric u-d QM to d quarks, for given p and ps. The other 
one can be analytically deduced from the following definition, 
i.e., 
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Here, the Fermi momenta of u and d quarks are expressed, 
respectively, as 
vu = (1 — 6/3)3y, 
va = (1+ 6/3)3y, (19) 
where v is the quark Fermi momentum of the symmetric u-d 
QM at the quark number density. p = 2p,, = 2pq. 
According to Eq. (18), Eq. (17) can be explicitly written 
as 
Esym(Pv, ps) = 
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where m = mo + Dp, / +0rn pl with Mmo = Muo = Mao. 
Here, we have assumed Myo = Mao. If Muo Æ Mao, then 
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We found that the discrepancy between the symmetry ener- 
gies given by Eq. (20) and Eq. (22) is essentially negligible. 
Therefore, in the following calculations, we adopted Eq. (20) 
to calculate the SQM symmetry energy. 
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Fig. 1. (Color online) Stability windows for the SQM, for different values of the isospin strength parameter C1. The absolute-stability region 
between the black solid and dashed lines widens significantly with increasing Cj, yielding results consistent with both the astrophysical 
observations of massive compact star PSR J1614-2230 with 1.928 + 0.017 Mo and tidal deformability 70 < Ai.4 < 580, measured for 
GW170817. The blue line denotes the maximum mass of the SQSs with 2 Mo, the red line represents the upper limit of tidal deformability 
measured for GW170817, i.e., Ai.4 = 580, and the red line with circles in panel (c) gives the lower limit Aj.4 = 70. The typical model 


parameters are indicated in panels. 


B. Equations for calculating structures and tidal 
deformability of SQSs 


To calculate the structures and tidal deformability of SQSs, 
we briefly introduce relevant equations. To obtain the 
SQS structures, one should numerically solve the Tolman- 
Oppenheimer- Volkov (TOV) equation: 


dP GmE (1+ P/E) (1 + 4rr3P/m) 


= 2 
dr r? 1 — 2Gm/r ? (25) 
with subsidiary conditions: 
d 
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where G = 6.707 x 10745MeV7~? denotes the gravitational 
constant. 


Dimensionless tidal deformability is defined as [58-61] 
(27) 


where C = UG is the compactness and kə is the gravita- 
tional Love number with its value usually in the range of 0.2— 
0.3. To calculate A as a function of the star mass or radius 


R, we need the value of k2, which is defined as 
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where yr = y(R) and y(r) is obtained as the solution to the 
following differential equation: 
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Here, F (r) and Q(r) are, respectively, defined as 
_ 1-4nr?le(r) — P(r)|G 
Hiz a (30) 
and 
_ An e(r) + P(r) 
Q(r) = Fir) seca +9P(r)G + ey 
6 m(r) tAr? P(r) o i ea 
m eyed 


where v? = OP(r)/0e(r) denotes the squared sound veloc- 
ity. In the case of a bare SQS, the energy density on the star 


surface is finite, which requires a correction to be added to 
the calculation of yg to account for the energy discontinuity, 
reading 


(32) 


where e, is the energy density difference between internal and 
external regions. 


I. NUMERICAL RESULTS AND DISCUSSION 


To study the isospin effects on the EOS of the SQM and the 
SQS structures, the values of the model parameters should be 
designated. First, the current masses of quarks are given as 
Muo = 5 MeV, mao = 10 MeV, and mso = 100 MeV. If Eq. 
(20) is used for computing the symmetry energy, the values of 
the up and down quark current masses are Muo = Mao = 7.5 
MeV. The values of the parameters C, D, and C; are deter- 
mined by studying the SQM stability window, combined with 
the astrophysical observation data of both massive compact 
star PSR J1614-2230 with 1.97+0.04 Mo and tidal deforma- 
bility 70 < Ay.4 < 580 measured for GW170817. 

Shown in Fig. | are the SQM stability windows for differ- 
ent Ct, with the value labeled in each panel. The black solid 
line in each panel indicates 930 MeV per baryon for the SQM; 
the SQM becomes unstable if the energy per baryon exceeds 
930 MeV. Similarly, the dashed line in each panel corresponds 
to the 930 MeV (per baryon) level for the u-d QM; below this 
level, the energy per baryon of the bulk u-d QM is smaller 
than 930 MeV, and nuclei are likely to decay to the u-d QM. 
This case contradicts the modern nuclear physics and should 
be forbidden. Therefore, only the area between the black 
solid and dashed lines is the SQM’s absolute-stability region. 
In addition, the blue line in each panel represents the config- 
uration that can yield the maximum SQS mass with 2 Mo. 
However, when Cy and/or C are negative, and baryon number 
density is sufficiently large, further decreasing C1 will lead to 
a negative quark mass and the requirement of maximum mass 
with 2 Mo will be no longer satisfied. Therefore, the blue 
lines in Fig. | will end at a certain point that is indicated by 
a star. The red line indicates the upper limit of Ay.4 = 580 
measured for GW170817. Moreover, it is worth mentioning 
that, if we decrease the values of the maximum mass Mmax 
or tidal deformability Aj.4, the lines representing them will 
shift upward, as indicated in panel (a) with up arrows. 

Comparing the three panels in Fig. 1, it is evident that, by 
increasing the isospin dependence parameter, C{, the dashed 
line will move significantly downward, whereas the solid line 
will be barely affected. Consequently, this will widen the 
SQS stable region. However, comparing panels (b) and (c), it 
can be observed that, for sufficiently high isospin parameter 
Ci values (e.g., for Cy = 400 MeV-fm°), further increasing 
it does not affect the absolute-stability region. In particular, 
we verified that beyond Ct > 10000 MeV-fm°, the absolute- 
stability region barely changed as the value of the parameter 
was further increased. 


By contrast, the blue lines for massive SQSs with 2 Mo 
and red lines for the SQS tidal deformability with A;,4 = 580 
tend slightly downward as C{ increases, but neither is affected 
significantly by varying Cj. To satisfy these two astrophys- 
ical observations (i.e., Mmax = 2 Mo and A14 = 580), 
the cross-points of the blue and red lines should be located 
in the absolute-stability region. Unfortunately, from these 
three panels, we can easily observe that, although the cross- 
point approaches the stable region with increasing C1, it does 
not end up in the stable region, even for Cy values as high 
as 10000 MeV-fm?. Given that the absolute-stability region 
does not change significantly by increasing Cy when it is al- 
ready sufficiently large, we conclude that in the equivparticle 
model with mass scaling as given by Eq. (2), the two observ- 
ables Mmax = 2 Mo and Aj.4 = 580 can not be simultane- 
ously satisfied by one set of model parameters. 


However, if we loosen the constraint on tidal deformabil- 
ity to its lower limit, namely decreasing Aj.4 to 70, the red 
line moves upward to the position of the red line with cir- 
cles, as shown in panel (c). Therefore, it can be inferred 
that as A4 4 decreases, the intersection point of the blue and 
red lines moves upward along the blue line. Therefore, we 
can choose the point (C, D!/?/MeV) = (—0.6, 164.5) in- 
dicated by the black dot in panel (c) as one of the typi- 
cal parameter sets, which is located in the SQM absolute- 
stability region. For comparison, we select parameter sets 
(C, D'/?/MeV) = (—0.7,168) and (0.08, 139) as typical 
parameter sets for investigating the SQM properties and the 
SQS structures, which are respectively indicated with a red 
square and blue diamond in panel (c). To study the effect of 
isospins on the SQM properties and SQS structures, we also 
chose points (C, D!/?/MeV) = (—0.6, 164.5) in panels (a) 
and (b) as typical model parameter sets indicated by black 
dots, which are in the SQM unstable region. 


In Fig. 2, the density behavior of the energy per baryon 
is shown for typical parameters. The figure reveals that 
in the equivparticle model, the minimum energy (circles) 
for each curve is at exactly the same density as zero pres- 
sure (asterisks), satisfying the requirement of thermodynamic 
consistency. Moreover, the one-gluon exchange interaction 
(C < 0) can significantly lower the minimum energy per 
baryon, which makes SQM much more stable than normal 
nuclear matter. In addition, comparing the three black lines, 
we find that with increasing Cj, the minimum energy per 
baryon increases. However, for high values of Cj, the dif- 
ference between the solid black line (C; = 10000 MeV-fm?) 
and dashed black line (Cj = 400 MeV: fm?) becomes small. 
This is because, with increasing C1, the difference between 
u and d quark fractions becomes small in -equilibrated 
SQM, which results in vanishing isospin asymmetry 6. We 
can examine this point in Fig. 3. With the parameter set 
(C,D'/?/MeV) = (—0.6, 164.5), the quark fractions are 
expressed as different Cy values indicated in each panel in 
Fig. 3. Evidently, the differences between u, d and s quark 
fractions become smaller as Cj increases. In particular, for 
Cı = 10000 MeV-fm~°, the difference between the u and 
d quark fractions is small, which leads to vanishing ô. Ad- 
ditionally, in this case, it is interesting to note that when 
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Fig. 2. (Color online) Density behavior of SQM energy per baryon 
in 8-equilibrium. It is obvious that the minimal energy (circles) 
for each curve is exactly at the density corresponding to zero pres- 
sure (asterisks), satisfying the requirement of thermodynamic con- 
sistency. 


P» = 0.6 fm~%, the u, d, and s quark fractions are essen- 
tially equal and approach 0.33, similar to the case of the color- 
flavor-locked (CFL) phase. The results for the density behav- 
ior of quark fractions in this study resemble those reported in 
Ref. [56]. 

Figure 4 shows the density behavior of the symmetry en- 
ergy for typical parameters. From the left panel of Fig. 4, 
it can be seen that the value of C plays an important role 
in determining the symmetry energy. Generally, for the nu- 
clear matter symmetry energy there are essentially two groups 
of density behavior, where one increases monotonously with 
density, and the other first increases and then decreases with 
increasing density [57, 62]. The latter is typically regarded as 
soft. Specifically, if the symmetry energy of nuclear matter 
becomes negative at suprasaturation densities, neutron stars 
exhibit pure neutron matter cores [63]. Similar to nuclear 
matter symmetry energy, the behavior of the SQM symmetry 
energy in the present work can be divided into two groups, 
depending on the value of the model parameter C. For non- 
negative C, the symmetry energy increases with increasing 
density. Otherwise, it first increases and then decreases with 
increasing density, which means that it becomes soft at high 
densities. However, unlike the nuclear matter symmetry en- 
ergy, the quark matter symmetry energy in our model does 
not become negative, because the quark masses given by Eq. 
(2) become negative at sufficiently high densities with nega- 
tive C. However, the soft quark matter symmetry energy will 
affect the chemical composition of the SQS core, and may 


Table 1. With the typical parameters (C, D'/?/MeV), the maxi- 
mum mass Max, the radius corresponding to the maximum mass 
R, the central density no, and the quark matter symmetry energy at 
the SQS center Esym are listed in the following table. 


(C, DY? MeV) |(0.08,139) | (-0.6,164.5) | (-0.7,168) 
Mmnax/Mo 1.93 1.97 2.01 
R (km) 11.39 10.44 10.61 
no (im *) 0.98 1.20 L11 
Esym (MeV) | 7493.97 | 874.38 | 215.93 


have important implications for the structures and dynamical 
evolution of SQSs; this issue requires further studies. 


Moreover, as will be discussed later, to satisfy the require- 
ment of the two-solar-mass constraint of the SQS maximum 
mass, the symmetry energies shown in the left panel can be as 
high as 3500 MeV, as indicated by the solid black line. In ad- 
dition, in the left panel of Fig. 4 we also show the results from 
Fig. 1 in Ref. [56]; these are shown here by the red dashed 
line indicated with DI-300x30. The reason for choosing the 
parameter set DI-300 for comparison was that the SQM sym- 
metry energy with DI-300 is of the same order of magnitude 
as the nuclear matter symmetry given by the RMF model with 
interaction NLpô. Therefore, from this panel, we can see that 
the SQM symmetry energy is approximately 30 times higher 
than the nuclear matter symmetry energy, which is consistent 
with the findings in Ref. [56]. In addition, it is interesting 
to note that at relatively low densities, the symmetry energies 
shown by the black solid and dashed black lines are indis- 
tinguishable. The right panel of Fig. 4 shows that, on one 
hand, when C; = 0, the symmetry energy increases with in- 
creasing density, while when Cy = 400 or 10000 MeV-fm?, 
the symmetry energy first increases and then decreases for 
densities above 0.4 fm~3; on the other hand, the symmetry 
energy increases drastically with increasing strength of the 
isospin dependence parameter Ct, e.g., the symmetry energy 
for Cr = 10000 MeV-fm? is approximately 20 times that for 
Cı = 400 MeV-fm®, or even more than 150 times that for 
vanishing Cy. 

Left and right panels of Fig. 5 show, respectively, the SQS 
structures and tidal deformability. The left panel of Fig. 5 
shows that, as the isospin dependence parameter C1 increases, 
the maximum mass and the corresponding radius of the SQS 
decrease, which is different from some previously reported 
results [56]. Moreover, although all maximum mass values in 
the left panel are above 1.93 Mo, namely, the lower limit 
of the observed mass of PSR J1614-2230, only the black 
lines with parametric sets are in the SQM absolute-stability 
region, which is illustrated in Fig. 1. This finding reflects 
the crucial role of the symmetry energy on the SQM EOS in 
the equivparticle model. Typically, the black solid line with 
(C, D'/?/MeV) = (—0.6, 164.5) and C; = 10000 MeV-fm® 
yields the maximum mass Mmax = 1.97 Mo, correspond- 
ing to the measured mass of PSR J1614-2230. In addition, 
owing to the fact that with increasing SQS central density the 
quark equivalent mass given by Eq. (2) becomes negative, 
the dashed black line where (C, D'/?/MeV) = (—0.7, 168) 
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Fig. 3. Quark fraction as a function of the baryon number density in the equivparticle model, with the corresponding parameters indicated in 
each panel. The blue, black, and red lines are respectively for the u, d, and s quark fractions. See text for a detailed discussion. 
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Fig. 4. Quark matter symmetry energy as a function of the baryon number density. The left panel shows that for C < 0 the symmetry energy 
peaks at a certain density. In addition, the red dashed line with DI-300 is replotted from Ref. [56], which has the same order of magnitude 
as nuclear matter symmetry energy. For the sake of comparison with the presently obtained results, we multiplied it by 30. Therefore, 
for satisfying the constraints imposed by astrophysical observations, the SQM symmetry energy can be as much as 30 times the nuclear 
matter symmetry energy. From the two panels we find that by increasing the isospin dependence parameter Ct, the symmetry energy can be 
drastically increased, especially for C > 0. 
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Fig. 5. Mass-radius relations of SQSs and the relation between the dimensionless tidal deformability and the SQS mass, for different model 
parameters. The solid black line for the parameter set (C, D!/?/MeV) = (—0.6, 164.5) and Cr = 10000 MeV-fm™? yields the maximum 
mass Mmax = 1.97 Mọo, corresponding to the measured mass of PSR J1614-2230. Although the dashed black line does not reach the 
maximum mass for the configuration with the model parameter set (C, D!/?/MeV) = (—0.7, 168) and Cr = 10000 MeV-fm~°, it yields 
the maximum mass above 2.01 Mo, which is in accordance with the observed mass of PSR J0348 + 0432. 
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Fig. 6. Maximal SQS mass, as a function of C1. Unlike the maximum SQS mass, which decreases monotonously with increasing Cr, the 
variation in the maximum mass of the u-d QS is complicated, depending on the sign of the model parameter C. 


does not reach the maximum mass. However, the mass of the 
most massive star given by the black dashed line can exceed 
2.01 Mo, which is in accordance with the measured mass 
of PSR J0348 + 0432 [28]. If we constrain the parameters 
to the SQM absolute-stability region in panel (c) of Fig. 1, 
we obtain that the SQS maximum mass can be as large as 
2.02 Mọ in the present model, which is still within the lim- 
its of the measured mass of PSR JO740 + 6620 with values 
of 2.08 + 0.07 Me [29]. Nevertheless, for GW190814 [34] 
it was shown that the binary merger possesses a secondary 
compact object with the mass as high as 2.50—2.67 Mọ at 
the 90% credibility level, which is far beyond the maximum 
mass of 2.02 Mo in the present model. If the secondary com- 
pact object is a neutron star or a quark star, this may imply 
that other inter-quark interactions should be taken into con- 
sideration, or a novel quark mass scaling should be proposed 
for the equivparticle model. Specifically, Table I lists the typ- 
ical quantities and symmetry energy for the SQS core. From 
Table 1, all the maximum masses shown are approximately 
2 Mo, and with a positively valued model parameter C', the 
SQS maximum mass yields a lower central density. However, 
it is worth mentioning that the quark matter symmetry energy 
at the SQS center for positive C is much higher than that for 
negative C’, although the central density is lower. 


The right panel in Fig. 5 shows the relation between the 
dimensionless tidal deformability and the SQS mass, for the 
equivparticle model with different parameter sets in Fig. 5. 
The horizontal blue dashed line in the right panel of Fig. 
5 corresponds to the SQS with 1.4 Mo, whereas the two 
blue solid lines indicate the range of the tidal deformability, 
that is, 70 < A < 580, based on the improved analysis of 
GW170817 using LIGO and Virgo collaboration [38]. Note 
that the tidal deformability increases rapidly as the SQS mass 
decreases, which means that lighter SQSs can be deformed 
easier than heavier SQSs. Meanwhile, it is evident that all 
intersections of the typical lines and the blue dashed line are 
located between the two solid blue lines. Given the maximum 
masses in the left panel, the typical configurations with the 
parameter sets we have chosen in the present model can sat- 
isfy both the astrophysical observations of the massive com- 
pact star PSR J1614-2230 with 1.97 + 0.04 Mo and the tidal 
deformability 70 < Ay.4 < 580 measured for GW170817. 


Shown in Fig. 6 is the dependence of the maximum mass 
on Ci, for the u-d QS (left panel) and SQS (right panel). 
Comparing the two panels, it is evident that, for the same 
model parameters, the maximum u-d QS mass is higher than 
that of the SQS. In addition, unlike the maximum SQS mass, 
which decreases monotonously with increasing Cj, the vari- 
ation in the maximum mass of the u-d QS is complex. In 
particular, for nonnegative C the u-d QS maximum mass first 
increases with Cy at Cy < 100 MeV-fm® and then decreases 
with Cy; by contrast, for negative C the maximum u-d QS 
mass first decreases with increasing Cr to approximately 400 
MeV-fm?, and then increases with increasing Cr. From these 
two panels, it is evident that for sufficiently large C; (e.g., for 
Cı = 10000 MeV-fm?), further increasing C; will not signif- 
icantly affect the maximum mass. Furthermore, Fig. 6 shows 
that for (C, V D /MeV) =(-0.6,164.5), the maximum u-d QS 


Table 2. With the typical parameters (C, D'/*/MeV), for the 
u-d QS maximum mass Mmax, the tidal deformability A,.4, and 
the energy per baryon E /n, are listed in the following table, where 
Ci = 10* MeV-fm?. 


(C, DY? /MeV) |(-0.6,164.5) | (-0.6,170) | (-0.6,175) | (-0.6,180) 
Mmnax/Mo ~2.30 | ~2.16 | ~2.04 | ~1.93 
ha 720 ~520 | ~380 | ~280 
E Jn (MeV) ~934 ~965 | ~993 | ~1021 


mass can be as large as 2.3 Mo, exceeding the observed 
2.08 Mo for JO740 + 6602. Nevertheless, in this case, the 
tidal deformability of the u-d QS yields Ay.4 œ 720, not in 
the 70 < Ay.4 < 580 range. This is illustrated in Table II, 
where the typical parameters (C, D'/?/MeV), the maximum 
mass Mmax, the tidal deformability A;,4, and the energy per 
baryon E /n, are shown for the u-d QS. Moreover, the data 
in the table suggest that, if the value of VD is increased to 
170, 175, or 180, the tidal deformability will be located in 
a reasonable range, but the maximum u-d QS mass will de- 
crease as well. Specifically, for (C, VD/MeV) =(-0.6,170), 
the energy per baryon for the u-d QM will exceed 930 MeV, 
implying an unstable state for the u-d QM, although the pa- 
rameter set will yield proper values of Mmax ~ 2.16Mọ and 
Ay.4 520, satisfying the constraints of astrophysical obser- 
vations. Therefore, for determining whether the u-d QS can 
simultaneously satisfy all of the above-mentioned constraints, 
further detailed investigations are needed. 


IV. SUMMARY 


In this study, we investigated the isospin effects on the 
SQM stability window by introducing an isospin-dependent 
term into the quark mass scaling to account for the quark- 
matter symmetry energy, and further studied the EOS of the 
isospin-asymmetric SQM. We found that by increasing the 
isospin dependence parameter C1, the SQM absolute-stability 
region becomes significantly wider, therefore enabling sup- 
port of the massive compact star PSR J 1614-2230 with 1.97 + 
0.04 Mo and the tidal deformability 70 < Ay.4 < 580 mea- 
sured for GW170817. As for quark fractions, we found that 
with increasing C1 the difference between u and d quark frac- 
tions became negligible, yielding vanishing isospin asymme- 
try 0; in addition, as the baryon number density increases, the 
u, d, and s quark fractions approach approximately 0.33, sim- 
ilar to the CFL phase results. The symmetry energy properties 
were then examined by considering different values of C{. It 
is interesting to note that for C < 0, with increasing C1 the 
SQM symmetry energy first increases and then decreases, im- 
plying soft symmetry energy. In addition, the symmetry en- 
ergy increases drastically with increasing isospin dependence 
parameter Cj, especially for C > 0. Owing to the widen- 
ing of the SQM absolute-stability region by incorporating the 
isospin effects into the SQM EOS, the measured mass of PSR 
J1614-2230 with 1.97 Mọ and the dimensionless tidal de- 
formability 70 < Ay.4 < 580 measured for GW170817 could 


be well explained as SQSs by choosing proper model param- 
eter values in the SQM absolute-stability region. However, 
considering that the measured mass of PSR J0740 + 6620 can 
be as large as 2.08 Mo, the present equivparticle model still 
cannot explain more massive SQSs with masses exceeding 
2.02 Mo. This may imply that other SQM effects should be 
considered, e.g., the presence of u-d quark Cooper pairs (2SC 
phase) [51], or a novel quark mass scaling framework should 
be developed [56]. Finally, unlike the SQS maximum mass, 
which decreases monotonously with increasing C1, the varia- 
tion in the u-d QS maximum mass is uncertain, depending on 
the sign of the model parameter C. 
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